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In this paper the dynamics of the classical chiral QCDi currents is studied. We 
describe how the dynamics of the theory can be summarized in an equation of the Lax 
form, thereby demonstrating the existence of an infinite set of conserved quantities. Next, 
the r matrix of a fundamental Poisson relation is obtained and used to demonstrate that 
the conserved charges Poisson commute. An underlying diffeomorphism symmetry of the 
equations of motion which is not a symmetry of the action is used to provide a geometric 
interpretation for the case of gauge group SU(2). This enables us to show that the solu- 
tions to the classical equations of motion can be identified with a large class of curves, to 
demonstrate an auto-Backlund transformation and to demonstrate a non linear superpo- 
sition principle. A link between the spectral problem for QCD2 and the solution to the 
closed curve problem is also demonstrated. We then go on to provide a systematic inverse 
scattering treatment. This formalism is used to obtain the reflectionless single boundstate 
eigenvalue soliton solution. 
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1. Introduction. 



QCD2 with massless quarks in the fundamental representation of the gauge group 
has received considerable attention recently. This has occurred for essentially two reasons: 
the first is the claim using non-abelian bosonization that the theory should be integrable 
0] ; the second stems from the observation that for massless quarks in a representation 
with dimension larger than or equal to the dimension of the representation of a heavy 
probe charge, the theory exhibits screening. 

In 1 + 1 dimensions, integrability is usually associated with the existence of a large 
number of soliton-like solutions of the classical equations of motion which are recast in a 
Lax form, from which the existence of an infinite number of conservation laws follows. If 
these survive quantization, they normally imply an elastic, computable S-matrix with 
no particle production. 

There is, of course, considerable evidence of particle production [|J || in the large iV c 
expansion of QCD-i |J . It has been argued recently that this is not necessarily inconsistent 
with the existence of an integrable sector in the theory []7j . The issue of screening versus 
confinement for massless quarks in the fundamental representation of SU(N) has also been 
discussed recently [|] . 

In a previous communication || , we analyzed the classical equations of motion of the 
theory in the light cone gauge A- = and established that the once integrated equation of 
motion for the currents is of the Lax form, demonstrating the existence of a infinite number 
of conserved quantities. We discovered that the equations of motion possess a continuous 
symmetry (diffeomorphism invariance in the x~ coordinate) which is not a symmetry of 
the action @. Making use of this symmetry and specializing to SU(2), we established that 
all classical Frenet-Serret curves with a constant term in the torsion are solutions of the 
equations of motion, and constructed a single soliton solution. 

The original motivation of was to use these classical solutions to construct bilin- 
ears [|TT] relevant to semiclassical approximations of fermionic theories. This has been 
implemented for the Gross-Neveu model ||12|| and the Thirring model . 

In this article, we are only interested in exploiting further the rich mathematical 
structure that emerges in the study of the equations of motion of the classical currents of 
the model. Specifically, in the gauge A_ = these take the form 



1 Recently, a similar symmetry has been noted in II B supergravity|10| . 
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(1.1) 



where 



d 2 _ 



0~ab — 2jab, 



(1.2) 



and jab is an 



SU(N) current!. 



In section 2, we show how equation ( p. . 1|) is arrived at and point out its Lax form. 
Following Fadeev, we construct the fundamental Poisson relation and recast the once 
integrated equation of motion as a zero curvature condition. In section 3, the chiral dif- 
feomorphism in the x~ coordinate is exhibited and for SU(2) it is used to map the most 
general solution of the equations of motion to appropriate Frenet-Serret curves. Using 
Darboux coordinates a second zero curvature condition is obtained. The Auto-Backlund 
transformation is identified and is used to obtain the single soliton solution of . This 
Auto-Backlund transformation is then used to construct a non-linear superposition princi- 
ple. In particular, we show how the full set of soliton solutions of the Sine-Gordon model 
are a subset of solutions generated by this non-linear superposition principle. For these 
solitons, we show further the the fermion number is quantized. The spectral problem of 
the zero curvature condition is found to be closely related to the spectral problem used to 
solve the closed curve problem. In section 4 our single soliton solution is shown to be the 
reflectionless single bound state eigenvalue inverse scattering solution. This construction 
is straightforwardly generalizable to SU(N). 

2. Some general observations. 

We work in light cone co-ordinated!, in the axial gauge A- = 0. In terms of the 
quark fields tpa = [V'-cuV'+a] which are in the fundamental representation of SU(N) 
(a = 1, 2, ...N), the Lagrangian density is given by 




(2.1) 



2 These equations of motion also arise in certain compactifications of M(atrix) theory[14] . 



3 Our conventions are = -4= (t ± x) 
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The equations of motion for both A + and ip+ do not involve derivatives with respect 
to x + and thus are interpreted as constraints. As usual, the constraints may be solved, 
leaving an effective theory for the quarks (x = (x + ,x~)): 



C = iV2i/jJ (x)d+i/j a (x) + y J dy x 



x ijjj (x)if; b (x)\x~ -y-\i) h ] (x + ,y-)i) a (x + ,y~) - — J dy~ x 1 " > 

x ip a ] (x)ip a (x)\x~ -y~\ip b \x + ,y~)ip b {x + ,y~). 

where we have let m —>■ and ip-a ~^ 4>a- The dynamics of the classical system now 
follows from this Lagrangian in the usual way 

id+if} a (x) = ^=a ab (x)i(j b (x) (2.3) 

where 



(J 



ab{x)= I dy \x -y \[tp a (y )ipt(y ) - j^5 ab tf; c (x + , V )ip c \x + ,y )]. (2.4) 



In this last expression, all fields are evaluated at the same point in "time" x + . In 
the above equations, the if) fields are normal commuting functions. Note that from the 
definition ( |2.4j ) we find 



d 2 _a ab (x) =2(Mx)^1(x) - ^5 ab Mx + ,y-)i>c t (x+,y-)y (2.5) 

It is useful to note that the classical equations of motion ( |2.3| ) and (|2.4j ) arise as the 
classical dynamics of the Hamiltonian 

H=~ J d 2 xTr(8 2 _a)a, (2.6) 
together with the (equal x + ) Poisson brackets 



{^ a ( x )^l(y)} = ^-S ab S(x-y). (2.7) 



An explicit realization for the Poisson brackets is 



if, / 5a 53 53 5a \ 

{a,/3} = -= / dx yl-rr -tt^ T7~7~T — i — • 2. 
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2.1. Lax Representation and Conservation Laws. 

Using the equations of motion, the time dependence of tpaipl is computed as 

• 2 

d+(lpa1pt) = -^PatpUcb - (TacMbl ( 2 - 9 ) 

which may be written in the form of a Lax equation 

. 2 

d+d 2 _a= i -^ r [d 2 _a,a]. (2.10) 
v2 

This last equation is nothing but the conservation of the probability 2-current. How- 
ever, due to the fact that the trace of the — component of the probability current vanishes, 
the trace of any power of the + component of the probability current is conserved (with 
respect to the "time" x + ). Since this quantity is local, it furnishes an infinite number of 
conserved quantities, indexed by the continous label x~ . Of course, this conservation law 
arises from the usual global U(l) symmetry. What is suprising however, is that in this 
case, it is not the total charge (i.e. the integrated charge density) that is conserved, but 
rather the charge density itself. Note the close similarity between (|2.10|) and the dynamical 
equation for the Heisenberg spin chain 

d + S = ig 2 [d 2 _S, S], (2.11) 

where S may be expanded in the basis of generators of SU{2). A close connection 
between the Heisenberg spin chain and SU[2) chiral QCD2 will emerge when we perform 
a systematic inverse scattering treatment in a later section. 

Consider next the time dependance of d-0~ a b(x). From the definition fl2.4|) 

d+d-a ab (x) = d+ J dy~e(x~ - y~)ip a (y~)ipl(y~) 

(2.12) 



W f 

= J dy ~ e<yX ~ ~ y~)[ d - a i a \^- 

Now, noticing that 



d-[d-a(x), cr(x)] = [d 2 La(x) 1 o~(x)] 
we find, after integrating by parts and dropping the boundary term 
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(2.13) 



• 2 

d+d-a ab (x) = -^[cLa(x),a(x)U (2.14) 

which is again of the Lax form. This Lax pair again gives an infinite number of 
conserved quatities contained in the conservation law 

d + Tr{d-a{x)Y = (2.15) 

The conserved charges discussed above are all gauge invariant quantities. Note that 
the equations ( |2.14| ) and ( |2.1U[ ) are invariant under diffeomorphisms of the form x~ — ► 



f(x~) || . The action however is not invariant under this symmetry. This diffeomorphism 
invariance will be used in the next section to provide a geometric interpretation of the 
problem. 

Finally we comment that it is also a simple task to show that 

d+Tr(a n (x)) =0 (2.16) 

which again represents an infinite number of conserved quantities. The conserved 
quantities ( 2.15 ) and ( |2.16 ) are both non-local functionals of the original field variables. 



For a discussion on non-local conservation laws for the Heisenberg magnet, see [IS| . 
2.2. The r matrix. 

The existence of an infinite sequence of conservation laws does not automatically 
guarantee integrability of the classical dynamics. One must prove in addition, that these 
conserved quantities Poisson commute. A very elegant way in which this is demonstrated 
is due to Faddeev fl6|] . The idea is to introduce the r matrix, defined by the fundamental 
Poisson relation (FPR) 



{U(x, A) ® U(y, (i)} = [r(A - fj), U(x, A) ® I + I ® U(y, X)]5(x - y), (2.17) 
where we have introduced the notation 



{A^B}„ 



k\h 



{An, Bkm}- 



(2.18) 



Usually the FPR is used to prove that certain conserved monodromies Poisson com- 
mute. The treatment for chiral QCD2 is much simpler. Considering U(x : A) = ■^d 2 _a{x), 
and taking traces of the above equation immediatly implies that 

{Tr(d 2 _a)(x),Tr(d 2 _a)(y)} = 0. (2.19) 

Thus, if we can find an r matrix such that (|2.171 ) holds, then we have demonstrated 
that our classically conserved quantities Poisson commute. For the case of the gauge group 
SU(2), it is not difficult to verify that d 2 _a(x) obeys exactly the same Poisson bracket 
relations as the Heisenberg spin chain variables a a S a (x). In addition, the r matrix for the 
Heisenberg spin chain, where U(x,\) = ^a a S a (x) is known. Thus, the r matrix defined 
by ( |2.17| ) is nothing but the r matrix of the Heisenberg spin chain! In a natural basis, the 
r matrix is given by r = P j (A — jj) Jl6j where 



10 
10 
10 
L0 1 

2.3. Zero Curvature Condition. 

We begin by noting that the integrability condition of the linear equations 



(2.20) 



d-a ab (x)(f) b (x) = a(j) a (x) 



(2.21) 



d + (j) a = - l -^=a ab (x)(p b (x 



reproduces ( 2.14j) . Now, shifting x 
theory to compute 5(f) = rjd-(f) we are led to 



(2.22) 

x~ + r\ and using first order perturbation 



d-4> a = -^=x + d-a ab (x)(p b . 



(2.23) 



Clearly, the pair of equations (|2.23|) and ( |2.22|) , which constitute the linear problem, 
also have the diffeomorphism symmetry. The integrability condition of (|2.23|) and ( |2.22j ) 
now provide a zero curvature representation 



[d+ + % a ' d - + % X+d - a] 







(2.24) 



of ( |2.14| ) . This zero curvature representation is well suited to a systematic inverse 
scattering treatment, which is provided in section 4. 
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3. The case of SU(2): Some Preliminary Comments. 

In this section we specialize our discussion to the case N = 2. 



3.1. The Geometry of the Lax Pair. 

Expanding the Lax equation ( |2.14| ) in the basis provided by the Pauli matrices (a 1 



as 



°ab{x) = ^f(x)a l ab (3.1) 



yields 



a 2 

d+d-f* = -^e'^d-Pf (3.2) 

After using the diffeomorphism symmetry to change to the new variable s(x~), such 
that d s pd s p = 1, we are able to interpret s as the arclength of the curve traced out by 
the three co-ordinates (/ , / , / 3 ) at any given time x + . Given the arclength and position 
of the curve it is a simple matter to compute the explicit expressions 

(ei)* = d 8 f (3.3) 



(e 2 Y = -d 2 s f (3.4) 

K 



(esY = -e^ k d s pd 2 j\ (3.5) 

K 

for the normal, binormal and osculating normal and 

« = VWW (3-6) 



T 2 
H 



l e ^ k d s fd 2 s fd 3 J k , (3.7) 



for the curvature and the torsion. The dynamics of the three unit normals ej is easily 
found, upon using (|3.2| ) 

d + (e n )i = -V2g 2 e^ k P(e n ) k (3.8) 
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In addition to this, the propagation of the curve in arclength is described by the 
Frenet-Serret equations which read 



d s ei = ne 2 



(3.9) 



d s e 2 = re 3 - nei 



(3.10) 



d s e 3 = -re 2 . 



(3.11) 



3.2. A Second Zero Curvature Condition. 

Now, transforming to Darboux co-ordinates 



z i = — — — (ei), + (e 2 )i + [e 3 )i = 1 



l-(ei)i 

we are led to the pair of Riccati equations 



(3.12) 



K 



where 



d-z l = -iTZi + -(zf + l) 



d+zi = -iazi + ^(6 + ic)zf + -(b-ic), 



(3.13) 
(3.14) 



a = -if l d s f\ 



(3.15) 



and 



b= — 

KT 



fdlf _ _ d 2 f j d 3 fjf i d 2 f i + K 2 fdsf 



K 



(3.16) 



1 



K 



(3.17) 



The Riccati equations derived above, and their relation to the Frenet frame will prove 
to be vital for the construction of auto-Backlund transformations. It is simple to verify 
that this is simply the coset space representation z = ^ of the following pair of linear 
eigenvalue problems 
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d a Pi~ 




- 1 K 

2 ' 2 




Pi 


d s p 2 _ 




K -i T 
L 2 2 - 




A. 



UnP 



d+Pi 



i 

L 2 



2 fl 



■|(c-«6) 



(c - i&) 



-a 





"A" 




.ft. 



(3.18) 
(3.19) 



The integrability condition for this pair of equations implies^ 

d+r = V2g 2 d+K = 0. 



(3.20) 



This last result is remarkable and unexpected. It implies that any curve may be 
mapped into a solution of QCD 2 . 

The connection between soliton solutions and curves has been discussed by a number 
of authors [I7| . In particular, Lamb [IS] shows how the single soliton solution of the 
Sine-Gordon model can be related to a curve of curvature k = asech(as) and torsion r = 
constant. This curve, upon identifying r = \^2g 2 x + , provides a soliton solution of the 
equation ( j3.2| ) . Explicitely, we have 



fi 

h 



2a 



2g 4 x+ 2 + a 2 
2a 

2g 4 x+ 2 + a 2 
2a 



sech(as) sin(V2g 2 sx + ) 
sech(as)cos(V2g 2 sx + ) 
tanh(as). 



(3.21) 



2g 4 x+ 2 + a 2 ' 

It is a straightforward excercise to explicitly verify that these fi do indeed satisfy the 
equation of motion (|3.2|) . We also state the solution with constant curvature k and torsion 



1 g 2 rr+ 



fl 

h 
h 



4 +2^ o 2 SiU (\ ^ + ^T X+2s ) 

g^x +z + 2k z V 2 



2k 



g 4 x + 2 + 2n 2 COS l K2 + Y X+2s) 



(3.22) 



2 4- 

g x^s 



f g 4 x+ 2 + 2k 2 



These restrictions on the curvature and the torsion can also be obtained directly as the 
integrability condition for the equations of motion ( |3.2|) and the Prenet-Serret equations B . 
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3.3. An Auto-Backlund Transformation. 

In this section, we will concentrate our attention on the class of constant torsion 
curves. Of course, this is only a small subset of all possible solutions to our equations of 
motion, but as we shall see, it includes a number of interesting solutions. 

The geometrical interpretation of a Backhand transformation is as a transformation 
that takes a given pseudospherical surface! to a new pseudospherical surface. Moreover, a 
Backlund transformation takes asymptotic lines to asymptotic lines. The fact that asymp- 
totic lines on a pseudospherical surface have constant torsion, is the first hint that it is 
possible to restrict the Backlund transformation to get a transformation that carries con- 
stant torsion curves to constant torsion curves. That this is indeed possible is the result of 
a theorem due to Calini and I vey [|K)f : Assume that the osculating normal^), binormal 
(e2) and tangent (ei) vectors to a constant torsion (r) curve, together with the curves 
curvature («) are known. The curve with Frenet frame given by 

e' x = ei + (1 — cos9)sin(3(cos(3e2 — sinflei) + sin9sin(3es, (3.23) 

e 2 = e<i — (1 — cos9)cos(3(cos(3e2 — sinflei) — sin6cos(3e^, (3.24) 

e' 3 = cos9es + sin9(cos/3e2 — sinftei (3.25) 
where (3 solves the equation 

^ = Csin(3 - k, (3.26) 
as 

and £an(|) = ^, has torion r' = r, and curvature 

k = k - 2Csin(3. (3.27) 

From the results of the previous section, this clearly implies that the above transfor- 
mation constitutes an auto-Backlund transformation for QCD<i- 

Let us illustrate this Backlund transformation with an example. First, lets find the 
Backlund transformation of the "trivial" solution of ( |3.9|) - ( |3.11| ) , i.e. the solution with 
k = and r = \^2g 2 x + . A suitable Frenet frame for this curve is given by 

5 i.e. a surface with constant negative Gauss curvature. 
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ei 



" " 




sin(y/2g 2 x + s) 
cos(^/2g 2 x + s) 




cos(y/2g 2 x + s) 
—sin(\/2g 2 x + s) 







, e 2 = 


, e 3 = 


(3.28) 


-1 















For k = 0, (|3.26| ) has the solution 



P = -2arctan(e Cs+a ) 



(3.29) 



which corresponds to a curve with curvature k' = —2Csinf3 = 2asech(Cs + a). The 
tangent to this Backlund transformed curve follows directly from ( |3.23j ) 



ei 



-2a 1 



A T + 2 



a 2 +2g 4 x 
-2a 2 
a 2 +2g 4 x+ 2 



sech(as)tanh(as)sin(\/2g 2 x + s) — 
sech(as)tanh(as)cos(y/2g 2 x + s) + 



2V2a- 2 



9 x 



2a 



a 2 +2g 4 x+ 2 
2y / 2ag 2 x+ 
a 2 +2g 4 x+ 2 



cos(\/2g 2 x + s) 
sin(y/2g 2 x + s) 



a 2 +2g 4 x+ 2 



sech (as) 



(3.30) 



This is nothing but the soliton solution ( |3.21j ) we found earlier. 

It is clear that in order to implement the Backlund transformation we need to solve the 
associated equation ( |3.26|) , which is difficult for non trivial curvatures. In the remainder of 
this section, we describe how to solve ( [3.26 ) for any curve whose Frenet frame is explicitly 
known and can be continued analytically as a function of r, with fixed curvature. 

First, setting y = tan^, (|3.26| ) becomes 



dy k .„ 9 . 



(3.31) 



Now, if we analytically continue the (real) torsion to r = —ia, then we find that we can 
identify y = —Z[ where Z\ is a Darboux coordinate ( |3.12[ ) satisfying the Ricatti equation 
( [3.13|) . This provides a direct way of obtaining (5. 

Lets illustrate this last point with an example. From the Frenet frame (|3.28| ) we 
construct the Darboux coordinate Z2 = cos(rx) + isin(rx). Evaluating this at r = — ia 
yeilds y = —e ax , which implies (3 = —2arctan(e ax ) in agreement with ( |3.29|) . 
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3.4- Fermion Number. 

Our classical solutions should correspond to solutions with a well defined fermion 
number. Thus, we should be able to identify an integer valued charge associated with each 
of the above classical solutions. In this section we show that this fermion number arises in 
a very natural way. 

First, note that from (2.3), it is clear that the quark field ip a is an eigenvector of d 2 _a 
with eigenvalue equal to the fermion number. Since d_ a is a two dimensional traceless 
Hermittian matrix, it has two real eigenvalues of opposite sign. The positive eigenvalue 
corresponds to the fermion number density. Thus, we may write 

V£V« = yj\Tr{dla)*. (3.32) 

In the curve language, this is simply the curvature. Thus, the total fermion number 
should be proportional to the integral over arclength of the curvature of the solution. 

Now, consider the set of solutions that arise as Backlund transformations of the solu- 
tion with k = . These solutions have curvature 



K = -2Csin(p), (3.33) 

where 

^ =Csin{(3). (3.34) 
as 

To obtain a soliton solution, we seek solutions which have the property ^ -> as 
s — > ±oo . From the above equation, this tells us that (3 — > rm as s — > ±oo for any integer 
n . Thus, the total fermion number is proportional to 



J cIsk(s) = - J ds2Csin((3) = -2 J dsj- = -2(3]^ = -2(m - n 2 )n. (3.35) 

Thus, the fermion number for these solutions are an integer times 2n . Now, consider 
the solution obtained by Backlund transformation of the above solution. In this case 
the curvature is given in (3.27) where the dynamics of (3 is now (3.26). The boundary 
conditions for f3 are unchanged, so that the above argument again leads to the conclusion 
that the total fermion number is an integer multiplied by 2n. Indeed, any solution which is 
obtained from a sequence of Backlund transformations of the k = solution, have a fermion 
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number equal to 2tt multiplied by an integer. The above arguments can be generalized to 
solutions obtained using the auto-Backlund transformation on a curve with an arbitrary 
initial curvature. In this case, the fermion number of the solutions obtained through the 
Backlund transformation will differ from the fermion number of the original curve by 27m 
where n is an arbitrary integer. There is however no restriction on the fermion number of 
the initial curve. Note that the results discussed in this section depend on the fact that 
x~ is identified with the arclength along the curve. 

3.5. Closed Curves. 

In this section, we consider the following differential-geometric problem: If a curve 
in R 3 is closed, then the curvature and torsion are necessarily periodic. However, given a 
periodic curvature and a periodic torsion, we do not in general obtain a closed in curve. In 
this section, we will show that the extra constraints that must be placed on the curvature 
and the torsion so as to obtain a closed curve, are naturally expressed in terms of the 
spectral problem for QCD^. More precisely, the requirement that a curve is periodic 
can be expressed as a constraint on the spectral curve of classical chiral SU{2) QCDi- 
It is interesting to study this curve, since spectral curves associated with other (1 + 1) 
dimensional models have been used to derive the Seiberg-Witten differential (and therefore 
also the exact mass-charge relationship of the BPS saturated states appearing in the low 
energy dynamics [BOf ). The observations made in this section are based on the recent 



proof given in fl2]J where Riemann surfaces corresponding to periodic curves are identified. 

As a start, consider the zero curvature condition given in ( |3.18| ) and ( |3.19| ) . Since 
the zero curvature condition is invariant under gauge transformations, we can equivalently 
consider the spectral problem! 



U^U = g- 1 Ug-g- 1 d s g, 
V ^V = g- 1 Vg-g- 1 d x+ g, 
where the connection U, V are defined in ( |3.18|) and (|3.19|) , and 



(3.36) 
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exp(^ dx'r(x')) 
exp(— | dx'r(x')) 



(3.37) 



6 In what follows, we have normalized the length of the closed curve (given by traversing the 
curve once) to 2tt. 
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Although this gauge transformation changes the connections U and V, it leaves the 
integrability condition unchanged. Concentrate on the problem 



L(X)F( S ,X) = 0=(J^-u)F{a,X), L(X) 



d_ 

ds 



-iiA 
hiq(s) \iX 



hiq(s) 



(3.38) 



The crucial elements of the above problem are that the magnitude of q(s) is related to 
the curvature of the curve, and that the derivative with respect to arclength of the phase of 
q(s) is related to the torsion. It is this property of QCD2 that makes it possible to consider 
the closed curve problem. The link just described, between a curve and a function q(s) is 
known as the Hasimoto map ||17|| . In what follows, we will need to make use of the shift 



operator, defined by the action 



/( a )-/( S + 27T) 



(3.39) 



Here, since we are looking for a closed curve, both the torsion and curvature are 
periodic. It thus makes to sense to think that our quarks will be moving (at fixed x + ) 
in a periodic potential. In analogy to the treatment of waves moving through a periodic 
potential (lattice) we define Bloch functions as the functions which are eigenfunctions of 
both the shift operator and L(X) 



4 1} (s : X) 



s,X)\ 



if. 



(2) 



4 2 \s,X) 



4 2) ( S ,A)J 



(3.40) 



L(A)^ (1) (s, A) = 0, V (1) (s + 2tt, A) = u^(X)^(s, A), 
L(A)^ (2) (s, A) = 0, tp {2 \s + 2tt, A) = u^(X)^ 2 \s, A). 

Since U(s, A) is traceless, the normalization of the Bloch functions does not change 
with arclength s, so that we must have 



^ (1) (A)w (1) *(A) = u^(X)J 2 >(X) = 1. (3.41) 

A point in the comple plane A belongs to the spectrum of L(X) if and only if \uj^ (A) | = 
1. Thus, we have obtained our first characetrization of our spectral curve. Now, we will 
use the concept of a transition matrix, which is introduced in the next section. The reader 
is asked to consult that section for details.@ For a generic complex A, the transition matrix 

7 We trust that this will not confuse the reader. 
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T(A) will have two eigenvalues u^(X) and u/ 2 )(A) and a pair of corresponding Bloch 
functions. In terms of these two functions, we introduce a hyperelliptic Riemann surface! 
Y such that u^(X) = uj(pi) and u/ 2 )(A) = u(ji2) where /Ui and \i2 are the pre-images of 
the point A under the projection Y — > C. Note that any hyperelliptic Riemann surface 
obviously has a natural holomorphic involution which amounts to transposing the two 
sheets. The function 

p(p) = ^-lnu(fi) (3.42) 

Z7TZ 

is called the quasimomentum function. In what follows, the quasimomentum differ- 
ential 

dp(fj.) = ^-dp (3.43) 

dp, 

will play an important role. 

The structure of the Riemann surface Y (called the Bloch variety) has been studied 
in detail in |22| . We will need the following definitions: A complex point A is called 
regular if u^\\) 7^ u/ 2 )(A), and irregular if uj^(X) = u/ 2 )(A). There are three types of 
irregular points (1) Branch points, (2) Non-removable double points and (3) Removable 
double points. An irregular point Ao is called a branch point if in going around this point 
we move from one sheet of Y to another (i.e. the monodromy around this point is non- 
trivial). If the monodromy around an irregular point is trivial, then this point is a double 
point. If the Bloch functions are equal at the double point, then it is non-removable; if the 
Bloch functions are not equal at the double point, then the double point is removable. 

One last piece of notation is needed now. Introduce the new SU(2) valued frame 
OiE % n , related to the Frenet frame by 

Ei = ei, 

E2 = cos(8)e2 — sin(9)e3, 

E3 = sin(9)e2 — cos(9)es, (3.44) 
9 = J K(s')ds'. 



A Riemann surface is called hyperelliptic if it is a two sheeted ramified covering of the 
Riemann sphere. 
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The SU(2) valued coordinate of the curve is now given by 



f(a)=f(0)+ / ds'E^s'). (3.45) 
Jo 

We are now ready to state the theorem that provides the restriction on the spectral 
curve: Let q(s) be a complex valued smooth periodic function of one variable s, q(s) 7^ 
for all s, q(s + 2tt) = q(s), Aq an arbitrary real number and T(s) the corresponding curve 
constructed from the curvature and torsion specified in q(s). Then, (1) the matrix Ei(s) 
is periodic with period 2n, if and only if Ao is a double point of the Bloch variety Y and 
(2) the function T(s) is periodic with period 2n if and only if Ao is a double point of Y 
and dp(fxi) = = dp((i2) where fi% and \ii are the preimages of Ao under the projection 
Y — > C. We are content with the statement of the theorem, and refer the reader to JHJ 
for a proof. 

To summarize, we have managed to show that the spectral problem in QCD2 can be 
used to define Bloch functions for periodic potentials. The requirement that these periodic 
potentials actually correspond to closed curves is then easily expressed as a condition on 
the quasimomentum differential. 

3.6. A Nonlinear Superposition Principle. 

The usual approach to massive QCD2 is to begin by performing a non-Abelian 
bosonization. Then, classical solutions of the bosonic action are found and a semi classical 
quantization may be carried out. The static classical solutions of these bosonic equations 
satisfy the Sine-Gordon equation. This allows the construction of multi-soliton solutions 
using the well known nonlinear superposition associated with the Sine-Gordon equation 
||23j| . In this section, we show that the nonlinear superposition principle also plays a role 
for the chiral theory. This explicitly demonstrates the existence of multi soliton solutions, 
which is another signal of the theories classical integrability. 

Begin by considering the spatial part of the Backlund transformation for the Sine- 
Gordon equation 

d ^~^ = Csin(p + <p), (3.46) 
where ip is the one soliton solution and hence solves 

^ = simp. (3.47) 

OS 

16 



From this last equation, it is clear that the single soliton solution corresponds to a 
curve with curvature —2 sirup. Now, noting that 



+ — = Csin((p + dp) + 2^f = Csin(ip + do) + 2Csin(ip) (3.48) 

(Jo (Jo 

we see that we may identify f3 = <p + do to obtain a Backlund transformation to a curve 
with curvature —2sirup — 2sin(3. Now, lets compute the three soliton solution <p from the 
two soliton solution dp 

® = Csin(p + (p) . (3.49) 

Rewriting this equation as 

g(g + 0) _ Csin(ip + <p) + 2^- = Csin{<p + dp) + 2(Csin(tp + <p) + sinUp)), (3.50) 

OS OS 

so that identifying (3 = <p + <p provides a Backlund transformation to a curve of 
curvature —2sin(ip + <p + dp) — 2sin(ip + dp) — 2sin(<p). By now, the generalization to the 
n soliton solution is obvious. 

From the nonlinear superposition for the Sine-Gordon equation, we obtain the alge- 
braic formula 



(3 = 2arctan 



n • r-2, r :; 
tan\ 



(2^)] + ^ + (3 - 51) 



-Ci - c 2 

for the Backlund transformation contructed from the n + 1 soliton solution of the 
Sine-Gordon equation. In this last formula, (p^ is the n soliton solution (of the Sine- 
Gordon equation) constructed from the n — 1 soliton solution (p n -i (of the Sine-Gordon 
equation) with parameter c±, §\ is the n soliton solution (of the Sine-Gordon equation) 
constructed from 4> n -i with parameter c 2 . Thus, the nonlinear superposition principle for 
the Sine-Gordon equation can be used to construct the multi soliton solutions of chiral 
QCD 2 . 



4. Dynamical Analysis using the Inverse Scattering Method. 

In this section, we present a systematic inverse scattering treatment of the classical 
equations of motion for the gauge group SU(2). We reproduce the single soliton solution 
described in the previous section, clearly demonstrating the validity of the analysis. How- 
ever, in contrast to the previous section, all results in this section are easily generalized to 
SU(N). 
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4-1. The direct scattering problem. 

In this section, we study the linear problem (2.22) derived in section 2.3. The eigen- 
value problem of the Lax operator reads 

d-<p a = \d-a ab <j) b . (4.1) 



From the time dependance of d-0~, (|2.14| ) , it is clear that requiring that d+d-a — > 



as \x~ | — > oo implies that the commutator of a and <9_<r must vanish as \x~ \ — > oo. There 
are two ways in which this commutator may vanish. First, the field d-0~ itself vanishes. 
Since we are working in terms of the arclength variable s, we have 

(ei)i(ei)< = \Tr{d.a 2 ) = 1, (4.2) 

so that cLcr^Oas ^ oo. The second possibility is that a becomes proportional 
to a single Pauli matrix. Without loss of generality (thanks to global color rotation invari- 
ance), we may take a to be proportional to a z . The requirement fl4.2p forces <9_cr — » a z as 
|a; _ | — > oo, which fixes the relevant boundary condition for ( |4.1|) . This linear problem has 
been studied by Takhtajan |24| (see also Fogedby |25j ) in connection with the classical 
Heisenberg spin chain. In this section we will review the results of his work relevant for 
the present analysis. 

The two Jost solutions F(x, A) and G(x, A) defined by the boundary conditions 

F(x, A) — > exp(-iXasx) (4-3) 

as x — > oo and 

G(x, A) -> exp(-i\cr 3 x) (4.4) 
as x — > — oo, have the integral representations 

/•OO 

F(x, A) = exp(— iXasx) + A / K(x,y)exp(—i\a3y)dy (4-5) 

J X 

r 

G(x, A) = exp(— iXa^x) + A / N(x 1 y)exp(—i\o~3y)dy. (4-6) 

J — oo 

The kernels K(x,y) and N(x,y) are independant of the eigenvalue A and solve the 
Goursat problem 
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9K(x,y) dK(x,y) . . 

-a 3 +d-a(x) = (x < y), 



dx 



dy 



dx 

with the boundary conditions 



dN(x,y) , dN(x,y) 

-<j3+0-(j(x) - = {x > y), 



dy 



(4.7) 



(4.f 



d-a(x) — (T3 + iK(x, x) — id-cr(x)K(x, x)a 3 = 0, 



(4.9) 



d-a(x) — as + zN(x, x) — id-cr(x)N(x, x)a 3 = 0. 
These last two boundary conditions imply that 



(4.10) 



d-a(x) = (iK(x, x) — a 3 )a 3 (iK{x, x) — 03) , 



and 



(4.11) 



d-cr(x) = (iN(x, x) - cr 3 )a 3 (iN(x, x) - cr 3 ) _1 . (4.12) 

The linear system ( [4.1|) has the important property that if ^1 and \&2 are solutions, 
then ^2 = *&iA, with A a constant-^ This allows us to write 



G(x,X) =F(x,\)T{\), 



(4.13) 



where T(A), is the transition matrix. If ^ is a solution of (|4.1|) , then writing det(^f) 
as exp(Trln^) we find, upon using Tr(d-a) = 0, that det($) is independant of x~ . Also, 
since the Pauli matrices have the property that a* = —(?2&i&2, it is clear that if ^ is a 
solution of (|4.1| ) for a spectral parameter A, then 02^*01 is a solution for A*. Using these 
properties of the solutions to (|4.1| ) , it is not difficult to show that the transition matrix 
has the form 



T(A) 



a(A) -(b(X*Y 
6(A) (a(A*)* 



(4.14) 



where 



9 This property is easily proved by differentiating \& 2 with respect to x , and using ( |4.l| ) . 
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a(X)a(X*)* + b(X)b(X*)* = 1. 



(4.15) 



The scattering data for the operator d-O is the set s = {a(X),b(X);£j,rrij, Im(£j > 
0,j = l,...,n}. a(A) can be analytically continued to the half plane Im(X) > 0, £j are 
the zeroes of a (A), and rrij = d ^^Tj ■ This last expression is written with the implicit 
assumption that we consider only simple zeroes 

^.2. Time dependance of the scattering data. 

The transition matrix is independant of x~ . This allows us to extract the x + depen- 

2 

dance of the transition matrix, in the limit x~ — > oo. Now, recalling A = + f(x~), 

and using the boundary condition of the a field, we find 

l -^d+G = x~a 3 G. (4.16) 

Alternately, we may write G = FT, and use the boundary condition for F, which 
implies 

l -^d+G = x-a 3 e- tXa * x ~ + e - tXasX ~ ^d+T. (4.17) 

g 2 g 2 

Comparing these last two equations, we find that 

^ + (4.18) 



dx + dx + dx + dX 

This last equation has the solution 

T = T(A --£*+). (4.19) 

This is a little unusual. For most models (for example Sine Gordon model, Heisenberg 
spin chain, nonlinear Schrodinger model) the spectral parameter is a constant and the 
scattering evolve in time. Here we have a linear time dependance for the spectral parameter 
and a set of scattering data which do not evolve in time! 
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4-3. The inverse scattering problem. 

Takhtajan has shown the kernel K(x, y; t) may be reconstructed from the time depen- 
dant scattering data, using the Gelfand-Levitan-Marchenko equation0 



/"OO 

K{x~ ,y~;x + ) + $i(x~ + y~;x + ) + / K(x~, z~; x + )$ 2 (z~ + V~] x + )dz~ = (4.20) 
where x~ < y~ , 



-A* 
A 



(4.21) 



$ 2 = -i 



o a_A* 

d-A 



(4.22) 



and 



A(x ;x + ) 



1 



6(A ' X+) e^-d\ + f2 mjiX+) e* 

3 = 1 ^ 



2tv 7-oo Aa(A, x+ 



(4.23) 



Having solved (|4.20|) for a set of scattering data, d-O is then reconstructed using 

o • 

Soliton solutions: Reflectionless Potentials. 

If the reflection coefficient r(A) = is zero, then the kernel and inhomogeneous 
terms in the Gelfand-Levitan-Marchenko equation are reduced to finite sums over the 
discrete spectrum. In particular, 



N 



A(x ;x + ) 



3=1 



£.3 



(4.24) 



Now, making the ansatz 



K = M(x~,x~ 



o 



(4.25) 



10 Fogedby[25] has given a clear derivation of this equation in the context of the classical 
Heisenberg spin chain. 
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we find that Q4.20| ) reduces to a system of iinear aigebraic equations 



M(x~,x + ] 



(4.26) 



+ 



!H± P -ia(x~+y-) 



E 3 if e ~ J 





0. 



which may be solved explicitly In this way, we obtain the iV soliton solution for 
QCD 2 . 

For the one soliton solution, we take a single bound state eigenvalue £. In this case, 
the kernel which solves the GLM equation is given by 



K(x , z~ 





m 


2 (£-n 2 


\a 


2 |£| 2 [(£-£*) 2 e- lx ~ (? - e) - 


m 


2 e -ia>- (£-£*)] 



|a| 2 e r -i(x~£*-£ Z - €H 2 r -i((x-+£*z~) 
i(£.-£*) e m e 



X 



m* 



Inserting this kernel into ( |4. 1 1| ) yields 

(2/m(0) 2 



d-h = 1 



(2/m(0) 2 + (2ite(0) : 



i(e-e) 



m 



and 



where 



cL/x + id_h = y/l - (<9_/ 3 ) 2 ^ 



(4.27) 



(4.28) 



(4.29) 



arg{m) + 2Re(^)x + arctan 



Jm(0 



ki viei 2 -(Me)) 2 

tan/if2x _ /m(0 + ^| 2/ "' (s) 



(4.30) 



V |m| 
Now, from the time dependance of the spectral parameter, it is clear that m is a 
complex constant, and that 



£ = —=x^ + x^ + i-, 
y/2 2 



(4.31) 



where x + and a are real constants. Upon inserting this into ( |4.28| ) and ( |4.29| ) , we 
find that we are reproducing the one soliton solution discussed in section 3.2. 
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